Effect of Poisson ratio on cellular structure formation 
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Mechanically active cells in soft media act as force dipoles. The resulting elastic interactions are 
long-ranged and favor the formation of strings. We show analytically that due to screening, the 
effective interaction between strings decays exponentially, with a decay length determined only by 
geometry. Both for disordered and ordered arrangements of cells, we predict novel phase transitions 
from paraelastic to ferroelastic and anti-ferroelastic phases as a function of Poisson ratio. 



Predicting structure formation and phase behaviour 
from the microscopic interaction laws is a formidable task 
in statistical mechanics, especially if the interaction laws 
are long-ranged or anisotropic. In biological systems, the 
situation is further complicated because interacting com- 
ponents are active in the sense that informed by internal 
instructions (e.g. genetic programmes for cells) and fu- 
eled by energy reservoirs (e.g. ATP), they react to input 
signals in a complicated way, which usually does not fol- 
low from an energy functional. Therefore these systems 
are often described by stochastic equations Q, • One 
drawback of this approach is that typically the stochastic 
equations have to be analyzed by numerical rather than 
analytical methods. However, for specific systems struc- 
ture formation of active particles can be predicted from 
extremum principles. In this case, analytical progress 
might become feasible again, in particular if analogies 
exist to classical systems of passive particles. One exam- 
ple of this kind might be hydrodynamic interactions of 
active particles like swimming bacteria 3]. Here this is 
demonstrated for another example, namely mechanically 
active cells interacting through their elastic environment 



Our starting point is the observation that generation 
and propagation of elastic fields for active particles pro- 
ceed in a similar way as they do for passive particles like 
defects in a host crystal, e.g. hydrogen in metal 0, Q. 
For a local force distribution in the absence of external 
fields, the overall force (monopole) applied to the elastic 
medium vanishes due to Newton's third law There- 
fore each particle is characterized in leading order of a 
multipolar expansion by a force dipole tensor P^. For 
many situations of interest, including cells in soft media, 
this force dipole will be anisotropic and can be written 
as Pij = PriiUj, where n is the unit vector describing 
particle orientation and P is the force dipolar moment. 
The perturbation of the surrounding medium resulting 
from a force dipole P' kl positioned at r' is described by 
the strain tensor Uij(f) = djOlGikir^r^P^, where sum- 
mation over repeated indices is implied and Gij is the 
Green function for the given geometry, boundary condi- 
tions and material properties of the sample. The strain 
Uij generated by one particle causes a reaction of another 



particle leading to elastic interactions. The essential dif- 
ference between active and passive particles is that cells 
respond to strain in an opposite way as do defects. For 
defects the interaction with the environment leads to a 
linearized potential V = —PijUij 6]. For example, an 
anisotropic defect attracting the atoms of its host lat- 
tice turns away from tensile strain, in this way enhanc- 
ing the displacement fields in the medium. In contrast, 
for mechanically active cells like fibroblasts, experimental 
observations suggest that they adopt positions and orien- 
tations in such a way as to effectively minimize the scalar 
quantity W = PijU^ = —V For tensile strain, this 

implies that contractile cells actively align with the ex- 
ternal field. Because they pull against it, they reduce dis- 
placement. For fibroblast-like cells, this behaviour might 
have evolved in the context of wound healing, when cell 
traction is required to close wounds. For a translation- 
ally invariant system, the effective interaction potential 
between elastically interacting cells therefore follows as 
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where we have used d[Gik — —diGik- In general, W 
scales as 1/r 3 with distance r and has a non-trivial an- 
gular dependance. 

The simplest model for the elastic properties of the ex- 
tracellular environment is isotropic linear elasticity. Thus 
there are two elastic constants: the Young modulus E 
describes the rigidity of the material and the Poisson ra- 
tio v the relative importance of compression and shear. 
Its maximal value is v = 1/2 (incompressible material). 
If such a material is tensed in one direction, the shear 
mode dominates and it contracts in the perpendicular 
directions (Poisson effect). For common materials, the 
minimal value for the Poisson ratio is v = 0. Then the 
compression mode dominates and uniaxial tension does 
not translate into lateral contraction. Isotropic linear 
elasticity is a reasonable assumption for the synthetic 
polymer substrates commonly used to study mechanical 
effects in cell culture HtHl- I n the following, we there- 
fore use the Boussinesq Green function for particles ex- 
erting tangential forces on an elastic halfspace (that is 
translational invariance applies to two dimensions). An- 
alyzing Eq. (QJ then shows that W has a pronounced 
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FIG. 1: Elastic interactions of cells lead to string formation. 
Here we consider the elastic interaction of a string of aligned 
dipoles spaced at equal distance a with a single additional 
force dipole at horizontal distance x and vertical offset y, both 
for (a) parallel and (b) perpendicular orientation. 



minimum for aligned force dipoles for all possible values 
of the elastic constants 0, . Recently, such alignment 
of cells in soft media has indeed been observed experi- 
mentally [HQ. 

In the following we show that much physical insight 
into cellular structure formation can be gained by start- 
ing from the finding that elastic interactions of cells favor 
string formation. We first consider an infinitely extended 
string of aligned force dipoles spaced at equal distance a. 
An additional dipole is placed at a horizontal distance x 
and with a vertical offset y, compare Fig. To simplify 
our notation, we use non-dimensional quantities: energy 
W is rescaled with P 2 /Ea 3 and length with a. Since the 
medium is assumed to be linear, the superposition prin- 
ciple applies and the effective interaction potential can 
be found by summing up all pairwise interactions. For 
parallel orientation (Fig.UJi), we have 
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and for perpendicular orientation (Fig. QJ)), we have 
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Eqs. (0 and © can be analyzed further with methods 
from complex analysis. Briefly, the Poisson sum formula 
Y^=-oo f(^ n ) = f c dz coth(7r z)f(z) allows to turn 
the sums into integrals in the complex plane 13]. By 
bending the contour C around suitable branch cuts in- 
troduced by the singularities of the corresponding f(z) 
at z = ±x + iy, the integral forms of Eqs. @ and © can 
be evaluated in the limit of large x: 
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FIG. 2: Interaction VK" between a finite string of N aligned 
dipoles and another dipole as a function of distance x for 
N = 200, offset y — and Poisson ratio v — 0. The full line is 
the numerical evaluation of Eq. (0). The dashed lines are the 
different scaling laws derived in the main text. Inset: VK" for 
an infinite string for x < 1, y = and v — 0, 0.2, 0.3, 0.4, 0.5 
from top to bottom. 
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Thus, despite the long-ranged character of the elastic 
pair-interaction, the effective interaction between an in- 
finite string and a single dipole (and a second string, 
respectively) is short-ranged and decays to leading order 
~ ^Jxe~ 2i:x . There is one exception, namely v = 0, when 
W" ~ e -27rx J ^ £ Qr p ara ^ e i dipoles and W 1 - = for per- 
pendicular dipoles. In dimensional units, the exponential 
decay occurs on a length scale a/2ir set by the dipolar 
spacing a only, independent of the elastic constants. 

In practice, strings will be finite. Therefore we now 
consider a single dipole with parallel orientation inter- 
acting with a finite string of N dipoles aligned along the 
y-axis and centered around the origin. For simplicity, we 
further specify to offset y = and Poisson ratio v = 0. 
For N = 1, the two dipoles interact via W" = 1/(ttx 3 ). 
For TV = 200, the full interaction potential as obtained by 
numerical evaluation of Eq. is shown as solid line in 
Fig. [21 Clearly there are different scaling regimes, which 
can be explained in the following way. For x > N, the 
dipoles in the string cancel each other and the mechanical 
action of the string is equivalent to two opposing forces 
d=P placed at ±iV/2. Therefore the string effectively acts 
like one dipole of total strength NP and the asymptotic 
potential decays with the dipolar power law N/(ttx 3 ). 
For 1 < x < TV, the string can be assumed to be infinite, 
but the discrete spacing between the dipoles can still be 
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neglected. Then we can use the analogy to the electro- 
static problem of an infinite, homogeneously charged line. 
Using Gauss' law and the fact that the force density van- 
ishes gives that now the potential W" has to be constant. 
The height of this plateau is fixed by the boundary con- 
ditions, that is in our case by matching it to the large 
distance regime at x ~ TV; therefore W" ~ 1/N 2 . For 
x « 1, the finite spacing between the dipoles becomes rel- 
evant and the exponential decay 87re -27nc / \fx predicted 
by Eq. (@J) becomes valid. Finally, for x < 1 the interac- 
tion with the string is dominated by interactions with the 
closest dipole in the string and W" crosses over to the 
dipolar power law 1/(ttx 3 ). Fig. [21 demonstrates the nice 
agreement between this scaling analysis and the full nu- 
merical result. While these four scaling regimes are valid 
in general, the exact details are very sensitive to dipole 
orientation, off-set y and Poisson ratio v. In particular, 
variations in v can qualitatively alter the interactions, es- 
pecially for x < 1, as demonstrated in the inset of Fig. 
for an infinite string. In general, for an infinite string the 
large distance regime and the height of the plateau van- 
ish, and we are left with the exponential decay derived 
above. 

String formation is also known forpassive particles, 
most prominently for electric dipoles [lj, LUj . The short- 
ranged nature of the effective interaction between strings 
due to screening inside the string is well-known in this 
case [ljl, [HI [HI 5 but to our knowledge has not been dis- 
cussed before for force dipoles. It has several important 
implications for cellular structure formation. First, it 
suggests that long-ranged effects do not dominate struc- 
ture formation at particle densities sufficiently large as 
to allow formation of strings of aligned dipoles. Sec- 
ond, for high particle densities the strong dependence of 
the string-string interaction on v suggests that structural 
changes might occur as a function of Poisson ratio. 

In order to address these issues, we used Monte Carlo 
simulations to study cellular structure formation on elas- 
tic substrates as a function of reduced cell density p and 
Poisson ratio v. Because cells are supposed to interact 
only elastically, a circular exclusion zone was attributed 
to each cell and positional degrees of freedom were fixed 
at random, thus avoiding cell-cell contact (because we 
are in two dimensions, p < 0.907). We then relaxed the 
orientational degrees of freedom using a Gibbs ensemble 
on W with a small finite effective temperature, which can 
be interpreted as an stochastic element in cellular deci- 
sion making. Similar procedures have been used before 
for modeling cellular structure formation [l^. Based on 
our simulations we predict three different phases. At low 
density, we always find many short strings with hardly 
no correlation. Presumably due to the screening effects 
described above, no global order appears, except when 
an external field is applied. Therefore this phase might 
be termed paraelastic. At high density and small values 
of the Poisson ratio v, spontaneous polarization occurs. 
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FIG. 3: Phase diagram for posit ionally disordered cells. At 
low cell density, an orientationally disordered (paraelastic) 
phase (p) prevails. At high cell density, orientational order 
sets in, with a nematic string-like (ferroelastic) phase (f) at 
low values of Poisson ratio v and a isotropic ring- like (anti- 
ferroelastic) phase (af) at large values of v. 



This phase of aligned strings might be termed ferroelas- 
tic and is characterized by a non- vanishing nematic order 
parameter S. As v is increased at high cell density, the 
nematic order parameter vanishes again for v c > 0.32. 
Now the local structure is ring-like rather than string- 
like. This phase results from the Poisson effect and might 
be termed anti- ferroelastic. In Fig. 03 we show typi- 
cal configurations and the phase transition line to the 
ferroelastic phase as measured in simulations by setting 
S = 0A. 

The important role of the Poisson ratio for cellular 
structure formation can also be demonstrated for regu- 
lar arrangements of force dipoles. Fig. 2] shows six can- 
didate structures identified by Monte Carlo simulations 
for dipoles on square and hexagonal lattices. To iden- 
tify optimal structures we calculate the interaction per 
particle by evaluating the corresponding lattice sums. 
For this purpose, we decompose the structures into A- 
and B-sublattices of parallel strings, as indicated by grey 
and white colors. The main contribution to W origi- 
nates from interactions within a string and is given by 
— 2£(3)(1 + v\ where ((z) is the Riemann Zeta-function 
with £(3) ~ 1.29. Because the interaction with adja- 
cent strings decays exponentially as shown by Eq. (@J), 
the interactions with adjacent strings are dominated by 
interactions with the nearest neighbor strings and the 
lattice sums converge quickly. The plots in Fig. 0] show 
that on hexagonal lattices, a phase transition from a fer- 
roelastic (h-a,b) to an anti-ferroelastic structure (h-c) 
occurs at v c — 0.32. On square lattices, W decreases 
strongly with v for the anti-ferroelastic structure (s-c). 
However, here lattice geometry stabilizes the ferroelas- 
tic phase (s-a) over the whole range of ^, such that the 
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FIG. 4: Structure formation of cells positioned on square (s) 
and hexagonal (h) lattices. Elastic interactions result in a 
competition between ferroelastic (a,b) and anti-ferroelastic 
structures (c) as a function of Poisson ratio v and lattice ge- 
ometry. 



ferro-antiferroelastic transition does not occur. This find- 
ing also implies that on substrates with v « 0.4 cellular 
structures might be switched between ferro- and antifer- 
roelastic phases simply by varying the geometry of cell 
positioning. 

In summary, we have shown that elastic interactions of 
cells lead to surprising and non-trivial structure forma- 
tion as a function of cell density and Poisson ratio. At 
low densities, the composite material of cells and elas- 
tic medium (the tissue) behaves as a paraelastic mate- 
rial: no global ordering occurs except for an externally 
applied field. Only for high cell density global order ap- 
pears, but the details depend strongly on Poisson ratio v. 
For small values of v, a ferroelastic phase appears, that 
is a macroscopic anisotropic force-dipolar field builds up 
that contracts the medium unidirectionally. This situa- 
tion is reminiscent of certain pathological situations when 
wound contraction detoriates into uncontrolled skin con- 
traction (contracture). For large values of the Poisson 



effect leads to an anti-ferroelastic phase, which is macro- 
scopically isotropic, but which shows pronounced local 
ordering into ring-like structures. 

Elastic substrates appear to be the ideal experimen- 
tal system to test our predictions. However, the typical 
physiological environment of tissue cells are the hydrated 
polymer networks of the extracellular matrix and future 
work has to show how our results carry over to this situ- 
ation. In particular, more details of the cellular decision 
making process have to be included in this case. Al- 
though then isotropic linear elasticity theory is certainly 
not sufficient anymore, one might expect that small val- 
ues of the Poisson ratio v correspond to cells probing 
single fibers (no Poisson effect), while large values of v 
correspond to cells probing more macroscopic regions of 
extracellular space, which are incompressible due to hy- 
dration effects. 
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